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The origin of magnetic fields in white dwarfs remains a fundamental unresolved problem in 
stellar astrophysics. In particular, the very different fractions of strongly (= 1 MG) magnetic 
white dwarfs in evolutionarily linked populations of close white dwarf binary stars cannot 
be reproduced by any scenario suggested so far. Strongly magnetic white dwarfs are absent 
among detached white dwarf binary stars that are younger than approximately 1 Gyr. In 


contrast, in semi-detached cataclysmic variables in which the white dwarf accretes from a 


low-mass star companion, more than one third host a strongly magnetic white dwarf!. Here 
we present binary star evolutionary models that include the spin evolution of accreting white 
dwarfs and crystallization of their cores, as well as magnetic field interactions between both 
stars. We show that a crystallization- and rotation-driven dynamo similar to those working 
in planets and low-mass stars” can generate strong magnetic fields in the white dwarfs in cat- 
aclysmic variables which explains their large fraction among the observed population. When 
the magnetic field generated in the white dwarfs connects with that of the secondary stars, 
synchronization torques and reduced angular momentum loss cause the binary to detach for 
a relatively short period of time. The few known strongly magnetic white dwarfs in detached 


binaries, including AR Sco’, are in this detached phase. 


The vast majority of close binary stars containing at least one white dwarf form through common 
envelope evolution’. The emerging detached post common envelope binary stars evolve towards 
shorter orbital periods driven by angular momentum loss? and eventually become semi-detached 
cataclysmic variable stars (CVs)°. Despite this clear evolutionary link, the fraction of strongly 


magnetic white dwarfs differs drastically between both types of close white dwarf binaries. 


We know more than 160 CVs with a strongly (larger than 1 MG) magnetic white dwarf’. 
In most of these systems the white dwarf accretes from its Roche-lobe filling low-mass star com- 
panion along the magnetic field lines. Due to interactions of the magnetic fields of both stars, 
the rotation of the white dwarf is synchronised with the orbital motion of the secondary star, such 


systems are called polars. In a smaller but still important fraction of CVs, the magnetic field of the 


white dwarf disrupts the inner parts of the accretion disk but the magnetic fields of both stars do 
not connect. In these so-called intermediate polars the rotation of the white dwarf and the orbital 
motion are therefore not synchronised’. A recent study of a volume-limited sample has shown that 
over a third of all CVs contain a magnetic white dwarf, with polars and intermediate polars making 


up approximately 29 and 7 per cent of the total population!. 


The situation is very different among the detached systems that are believed to be the pro- 


genitors of CVs. Only 15 strongly magnetic white dwarfs in close detached binaries are currently 
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known” and they make up less than two per cent of the over one thousand known systems 
All but one of these detached magnetic white dwarf binaries have been identified due to a pecu- 
liar emission line!?, which turned out to be the third harmonic of a cyclotron fundamental emitted 
by a low-density plasma. This observational finding is convincingly interpreted as resulting from 
wind accretion onto a strongly magnetic white dwarf in a close but detached binary'*. Because 
angular momentum loss will eventually drive them into a semi-detached CV configuration, these 
systems have been termed pre-polars’. The only exception is AR Sco, which was discovered be- 
cause of the optical and radio pulses of its rapidly (1.97 min) rotating white dwarf*. AR Sco and 
most pre-polars are located in the orbital period range between 3 and 5 hours, their secondary stars 
are close to filling their Roche-lobes (R/R typically exceeding 80 per cent)!°, and the magnetic 
white dwarfs are cold with effective temperatures typically below 10,000 K which implies that they 
formed at least one Gyr ago. Among younger close detached white dwarf binaries, not a single 


strongly magnetic white dwarf has been found!>!°. This puzzling situation, a large fraction (ap- 


proximately 35 per cent) of strongly magnetic white dwarfs in CVs and the much smaller fraction 


of magnetic white dwarfs among the detached systems that are believed to be CV progenitors, was 
aptly expressed by Liebert’? who asked Where Are the Magnetic White Dwarfs with Detached, 
Non-degenerate Companions? This fundamental question of stellar evolution has remained with- 
out an answer for 15 years and the solution must link the evolution of white dwarf binaries and 


magnetic field generation in white dwarfs. 


Several hypotheses have been put forward to explain magnetic field generation in white 
dwarfs in the last decades. In their present form, the fossil field scenario'’, a dynamo operat- 
ing during common envelope evolution!®, or coalescing double white dwarfs!’ face one serious 
problem when considered as the main formation mechanism for close magnetic white dwarf bina- 
ries. In all these cases the magnetic fields are generated during the formation process of the white 
dwarf, i.e. according to the current versions of these theories also detached post common envelope 
binaries containing hot white dwarfs should host strong magnetic fields which is clearly not the 


case!®, 


As the only mechanism that depends on the age of the white dwarf, a dynamo similar to those 
operating in low-mass main-sequence stars and planets” has been suggested to operate if the white 
dwarf is rotating rapidly when the core of the white dwarf is crystallizing??. As a white dwarf 
with a Carbon/Oxygen core cools, the ions in the core begin to freeze in a lattice structure. During 
crystallization the chemical potential, temperature, and pressure of the liquid and solid phases 
must remain equal. As a consequence, in crystallizing white dwarfs the solid phase becomes 


richer in oxygen and sinks while the carbon excess mixes with the outer liquid envelope which 


is redistributed by Rayleigh—Taylor instabilities. This configuration is similar to that found in the 
core of the Earth, where the light element release associated with the inner core growth is a primary 


driver of the dynamo generating the Earth’s magnetic field”!. 


The predicted field strengths generated by this dynamo in white dwarfs depend entirely on 
the scaling law that is assumed. Applying the same scaling law that is used for planetary magnetic 
fields leads to field strengths in white dwarfs of up to approximately 1 MG”. However, taking into 
account a likely dependence of the dynamo efficiency on the magnetic Prandtl number and that 
in white dwarfs this number is orders of magnitudes larger than in planets, one could expect this 


dynamo to be able to generated fields much stronger than 1 MG (for more details see Methods). 


Assuming that the crystallization and rotation driven dynamo can generate strong magnetic 
fields we computed evolutionary tracks starting with non-magnetic post common envelope bina- 
ries. We incorporated several physical mechanisms in the stellar evolution code MESA”. The 
spin-up of the white dwarf due to the accretion of angular momentum”? was included as well as 
crystallization of the white dwarf**. We furthermore incorporated the reduction in angular momen- 
tum loss through magnetic braking when the magnetic fields of both stars connect”, a condition 
for synchronisation”®, and a prescription for angular momentum transfer from the spin of the white 


dwarf to the orbit during synchronisation (see Methods for more details). 


Figure | shows an exemplary evolutionary track resulting from our simulations and a schematic 
illustration of the different stages the system passes through. In our model the white dwarf is born 


without a strong magnetic field which is in agreement with the absence of strongly magnetic white 


dwarfs in young close detached post common envelope binaries (stage 1). The non-magnetic de- 
tached binary star then evolves slowly towards shorter orbital periods and the white dwarf cools. 
When the secondary fills its Roche-lobe the binary becomes a CV, accretion spins up the white 
dwarf” and the core may be crystallizing (stage 2). If this happens, both conditions for the crystal- 
lization and rotation driven dynamo are met and a strong magnetic field is generated (stage 3). If 
the magnetic field is strong enough to connect with the field of the secondary star the latter provides 
a synchronising torque on the white dwarf spin. The spin angular momentum transferred to the 
orbital motion due to the synchronising torque causes the secondary to detach from its Roche-lobe. 
This effect is supported by the fact that the connection of the two magnetospheres reduces angular 
momentum loss due to magnetic wind braking because part of the gas escaping from the secondary 
is accreted by the white dwarf along the connected field lines, instead of leaving the system*>. At 
the beginning of synchronisation the system becomes a radio pulsing white dwarf binary, such 
as AR Sco, in which a rapidly rotating white dwarf spins down on a relatively short time scale 
(stage 4) of approximately 1 — 10 Myr. When the white dwarf spin is synchronised with the orbital 
motion of the secondary star, the binary star has evolved into the pre-polar stage in which wind 
accretion generates the observed cyclotron emission (stage 5). Angular momentum loss through 
reduced magnetic braking and gravitational radiation brings the stars closer together until the sec- 
ondary again fills its Roche-lobe and finally a polar is formed (stage 6). The detached phase (stage 
3-5) occurs preferentially for orbital periods shorter than roughly 5 hr as the reduction of magnetic 
braking and the transfer of spin angular momentum to the orbit have less impact at longer periods. 


Detached phases occur less frequently at orbital periods shorter than approximately 3 hr because 


the conditions for magnetic field generation are typically met at early stages of CV evolution (see 


Methods for a detailed comparison of observations and model predictions). 


While a different mechanism such as common envelope merger'® and/or double white dwarf 
merger’? is probably required to explain most single white dwarfs with magnetic field strengths 
exceeding 1 MG (see Supplementary Material for details), the new evolutionary sequence predicted 
by our model is in excellent agreement with the observations of close binary stars containing 
strongly magnetic white dwarfs. As the conditions for the dynamo are relatively frequently met in 
CVs (see Supplementary Material for a discussion of the predicted relative numbers), our model 
successfully reproduces the large fraction of CVs containing magnetic white dwarfs. The few 
known detached binaries containing magnetic white dwarfs are not progenitors of CVs but instead 
systems that have been non-magnetic CVs in the past and where mass transfer stopped because the 
field was generated. This naturally explains why there is not one single hot and young strongly 
magnetic white dwarf in a detached binary. Instead our model predicts that all of them are older 
than 1 Gyr with crystallizing cores as is observed (Fig. 2) and also explains why all 15 currently 
known close detached binaries containing magnetic white dwarfs are close to Roche-lobe filling 
(Supplementary Figure 1). Furthermore, the observed clustering of pre-polars in the orbital period 
range of 3 — 5hr is well reproduced in our simulations as the detached phase typically occurs in 
this range. Finally, our evolutionary sequence predicts short phases of rapidly rotating and strongly 
magnetic white dwarfs in detached binaries thereby offering an explanation for the existence of the 
radio pulsating system AR Sco. We conclude that a dynamo-driven generation of magnetic fields 


in white dwarfs resolves several long standing problems at once. 
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Fig. 1. Spin period, Roche-lobe filling factor, and mass transfer rate as a function of orbital period. 
The slowly rotating white dwarf (stage 1) in a non-magnetic post common envelope binary is 
spun-up due to accretion (stage 2). In case the white dwarf is crystallizing a strong magnetic field 
is generated by a dynamo (stage 3) and as soon as it connects with the field of the secondary the 
system becomes detached. At first it resembles the radio pulsing white dwarf binary AR Sco (stage 
4). When the spin of the white dwarf and the orbital motion are synchronised (stage 5) the system 
becomes a pre-polar before mass transfer starts again (stage 6). When the secondary becomes 
fully convective, a second short detached phase (between ~ 2.07 and ~ 2.09 Gyr) occurs which 


represents a short version of the classical orbital period gap (hatched region). 
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Fig. 2 Observed white dwarf effective temperatures versus their masses for non-magnetic post 


)'2.27. pre-polars (red)!°, 


common envelope binaries from the Sloan Digital Sky Survey (light gray 
and the first radio pulsing white dwarf binary AR Sco (blue)*:”°. The errors on the observations 
correspond to 1g and arrows indicate upper limits. The lines represents the mass fraction of 
crystallized matter in the white dwarf interiors. Solid, dashed, dash-dotted and dotted lines show 
the temperature for each white dwarf mass corresponding to the onset of crystallization, 10 per cent 
of mass is crystallized, 50 per cent and 80 per cent, respectively”. As predicted by our scenario, 
magnetic systems host on average much older and more massive white dwarfs than their non- 
magnetic counterparts. The trend towards larger masses in magnetic systems is totally consistent 
with the observed white dwarf mass distribution in cataclysmic variables?”*° which fits with the 


model prediction that these systems have been CVs in the past. In addition, the white dwarfs in 


AR Sco and all pre-polars are consistent with having crystallizing cores. 


Methods 


The convective dynamo and scaling laws. The magnetic fields of Earth and Jupiter, as well 
as those of rapidly rotating low-mass stars, are generated by convection-driven dynamos. The 
basic ingredients for these dynamos to work are a strong density stratification, and an extended 
convection zone. To derive the magnetic field strength from fundamental properties of a given 
planet or star the use of scaling laws is required. In general it is neither clear which quantities need 
to be considered in these scaling laws nor whether a single scaling law exists for all planets and stars 
whose magnetic fields are generated by the convective dynamo. In the last decades, several scaling 
laws have been proposed, most of them based on the assumption of a balance between Coriolis 
force and Lorentz force*!. In recent years a scaling law derived from geo-dynamo models*” which 
assumes the energy flux to be the dominant parameter has become popular. This scaling law has 
been shown to successfully explain the field strengths of Earth and Jupiter as well as those of 


rapidly rotating fully convective stars”. 


In the first application of the convective dynamo mechanism to white dwarfs this energy flux 
scaling law was assumed??. The vast majority of white dwarfs consist of Carbon and Oxygen. 
As a white dwarf cools, the ions in the core begin to freeze in a lattice structure**. The phase 
diagram of the carbon—oxygen mixture is of the spindle form**:*>, therefore, as soon as the white 
dwarf has sufficiently cooled and starts to crystallize, the solid phase becomes richer in oxygen and 
sinks while the carbon excess mixes with the outer liquid envelope. As a consequence, an oxygen- 


rich core forms which is significantly denser than the carbon-rich liquid which is redistributed 
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by the Rayleigh-Taylor instability**?7 


. Such a strong density stratification in combination with 
convection represent the key ingredients for the dynamo generating the Earth’s magnetic fields and 
those of low-mass stars. Assuming the same scaling law that has been successfully applied to the 
Earth, Jupiter, and rapidly rotating fully convective stars for crystallizing white dwarfs predicts 
field strengths of at most a few MG”? which, if true, would exclude the dynamo from being the 


mechanism generating the much stronger 1 — 200 MG fields found in white dwarfs in close binary 


stars. 


Here we argue that this scaling law is most likely not applicable to white dwarfs. In the 
magneto hydrodynamic simulations that led to the energy flux scaling law, the magnetic Prandtl 
number (Pn, defined as the ratio between viscous and magnetic diffusivity) is assumed to be equal 
to one. However, in both fully convective stars as well as the Earth this number is significantly 
smaller (1074 — 10~° and 1076, respectively). To quote the discovery paper’? of the energy flux 
scaling law “we cannot rule out an additional dependence on other parameters, in particular the 
magnetic Prandtl number. Although the suggested dependence is weak, it poses a serious problem. 
Given the large range of extrapolation over five orders of magnitude from our models to planetary 
values of P», the results obtained from the scaling laws with or without a dependence on Pn 
differ substantially.” This early statement has been confirmed more recently by Brandenburg?’ 
who wrote: “This leads to a scaling law that has been verified over many orders of magnitude 
(Christensen et al. 2009). As we have seen, this scaling law must be affected by Pm...” . Given that 
the magnetic Prandtl number for crystallizing white dwarfs has been estimated to be 0.58%, i.e. 


roughly six orders of magnitude larger than that of the Earth and still 4 — 5 orders of magnitude 
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larger than that of fully convective low-mass stars, the dynamo could be more efficient in white 


dwarfs. 


In fact, the comparison between predicted and observed field strengths of single white dwarfs 
presented by Isern et al.” indeed indicates that observed fields in crystallizing single white dwarfs 
exceed those predicted by the energy flux scaling law. The energy flux scaling law is only appli- 
cable in the saturated regime of the dynamo, i.e. it only works for rapidly rotating objects with 
Rossby numbers below ~ 0.1. For white dwarfs, this implies spin periods shorter than ~ 90 s°. In 
white dwarfs that rotate substantially slower, the generated fields should be smaller than those pre- 
dicted by the energy flux scaling law. All the observed magnetic single white dwarfs within 20 pc 
(Table 1 of Isern et al.?°) rotate with periods ranging from hours to days. Despite rotating orders 
of magnitudes too slow for the energy flux scaling law to be applicable (which is applicable in the 
saturated regime), the field strengths of these white dwarfs even slightly exceed those predicted by 
the energy flux scaling law. This implies that, if a convection driven dynamo works in crystallizing 
and rotating white dwarfs, the energy flux scaling law underestimates the field strengths of slowly 
rotating single white dwarfs. The field strengths of rapidly rotating white dwarfs in the saturated 
regime should therefore as well largely exceed those predicted by the energy flux scaling law. A 
suitably revised scaling law that incorporates the dependence on the magnetic Prandtl number is 
therefore likely to predict fields much stronger than ~ 1 MG for white dwarfs that rotate in the 


saturated regime of the dynamo. 


Based on the above arguments, we here assumed that a rapidly rotating white dwarf (Pwp ~ 
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20 — 30s) with a fractional crystallized mass between 10 and 80 percent can drive a dynamo that 
generates the large magnetic field strengths observed in white dwarfs in cataclysmic variables (up 
to several 100 MG). We calculate the evolution of close detached white dwarf binaries and cata- 
clysmic variables (CVs) using MESA”, taking into account the spin-up due to accretion and crys- 
tallization due to cooling of the white dwarf as well as magnetic field generation if a crystallizing 


white dwarf rotates rapidly. 


Spin-up of the white dwarf. If a white dwarf accretes from an accretion disk in a CV its angu- 
lar momentum and therefore its angular velocity increases, while during nova eruptions the white 
dwarf looses mass and angular momentum. According to King et al.°, the resulting angular mo- 


mentum balance for an accreting white dwarf is given by 


dw . . 
r= = a(-M) (OMe Rae" +049 (Ahn Rim w), (1) 
where w is the white dwarf spin, J its moment of inertia, G the gravitational constant, M. > the mass 


transfer rate averaged over nova cycles, and Mwp and Rwp are the white dwarf mass and radius. 


The first term in the right-hand part of the equation corresponds to the spin-up due to accre- 
tion, while the second term represents the spin-down due to material leaving the white dwarf during 
nova eruptions. The parameter 0 < a < 1 represents the spin-up efficiency, which might be smaller 
than one due to mass loss processes not considered in equation Eq.[ij(e.g. winds from the inner ac- 
cretion disk or jets from the boundary layer/white dwarf). The parameter 0 < € = Mwp / Mp <1 
relates the mass loss of the white dwarf to the mass transfer rate. The range of possible val- 
ues reflects the fact that the white dwarf mass might decrease over a nova cycle, i.e. more mass 
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might be expelled during a nova eruption than was accreted between two eruptions. The parameter 
0 < 7 < 1 represents the square of the radius of gyration of the ejected envelope and depends on 


the geometry of nova eruptions. 


By solving the non-homogeneous differential equation [I] the WD spin as a function of the 
donor mass is given by 


C 
w= = + Cs exp (C2 Mo), (2) 
2 


where Mə is the donor mass in solar masses and the terms C1, C2 and C3 are given by 


E AÁ wko 
Ci =Q T2 Mens (3) 
C= — (146) +Be (4) 
k2 Mwp 
Ci 
C3 = Wo — Co exp (—Cy M20) (5) 
2 


where wpy is the break-up white dwarf spin, given by 


GMwp Mwp Rwa =| 
pe; 2 . 
Mn = an | (e) (Bm)? p, © 


The moment of inertia I is k? Mwp Rẹyp, with k given by 0.452 + 0.0853 logy) (1 — Mwp/1.44Mo), 


if Mwp < 1.368 Mo, and © 0.275, if Mwp Z 1.368 Mo. The constants A and B are given by 


3/4 
_ |1- 0.61 ( Mwno/Mo)*" @ 
1 — 0.61 (Mwp/Mo )*” 
p — _305 (Mwp/Mo pas] 0.051 (Mwp/Mo) (8) 


3 [1-0.61(Mwn/Mo)*]  ¥ [1 — 0-69 ( Mwn/Mo)] | 


The sub-index 0 in some quantities indicates their initial values just after the common-envelope 
evolution. 
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Synchronisation of the white dwarf. In accreting high-field magnetic white dwarf binaries the 
magnetic field dominates the accretion process. For sufficiently strong white dwarf magnetic mo- 
ments, the overflowing material is channelled along the white dwarf’s magnetic field lines as soon 
as the magnetic pressure dominates the gas ram pressure. For such strong white dwarf magnetic 
fields, the interaction of the magnetospheres of the white dwarf and the donor cause a synchroniz- 


ing torque on the white dwarf”. 


We assume here that the onset of synchronization occurs when the synchronizing torque is 
greater than the accretion torque (assuming that the magnetospheres of the white dwarf and the 


donor are connected)”®, i.e. 


aie > 111 (a) (ya) (m) (2) ">? ea [0.5 — 0.277 log, E) 

(9) 
where uwp is the white dwarf magnetic moment, Mə is the donor mass, Miot = Mwp + Mo is 
the binary total mass, P.r» is the orbital period and B» is the strength of the magnetic field of the 


secondary star. 


Once the condition provided by Equation [9] is satisfied, the spin-down process starts. For 


simplicity, we assume that the white dwarf spin decreases exponentially on a given time-scale, i.e. 
w(t) = Q + Cet (10) 


where 7 is the synchronisation time-scale, Q the orbital angular velocity, w the WD spin, and 


C = wo — 2 is a constant (because variations in Q are negligible during the spin-down phase). 
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This assumption gives 


dw (Q—w) 
dt ap ` H 


This equation of course represents a very crude approximation of a very complicated process. 
The torque on the primary as a function of asynchronism has been calculated previously in great 
detail“, but only for rather small asynchronism. Performing such detailed calculations of the 
magnetic field interaction for the large asynchronism in our systems, is beyond the scope of this 
paper. For the purpose of this paper, we assume that the total angular momentum loss by the white 
dwarf is transferred to the secondary, incorporated in the orbital motion, and that this occurs on the 


time-scale 7, 1.e. 


d (wI) = ez% 


dt T 


Jaa = >0. (12) 


In this simple model the torque is always decreasing with increasing synchronism which might not 


be correct? 


. However, this should not affect the validity of our model as only the amount of spin 
angular momentum transferred to the orbit and the time scale on which this happen are important 


for the secular evolution of a given system. 


Reduced magnetic braking. For secondary stars still containing a radiative core, angular momen- 
tum loss through magnetic wind braking drives mass transfer in CVs°!'. The strength of magnetic 
braking depends significantly on the extension of wind zones compared to dead zones. The latter 
correspond to regions where gas is prevented from escaping the star due to the pressure of closed 
magnetic field loops, i.e. the magnetic pressure is greater than the thermal pressure*!. If the white 
dwarf in a CV is strongly magnetic and the magnetospheres of both stars connect, the wind zone is 
42,43. 


reduced which leads to a decrease of orbital angular momentum loss due to magnetic braking 
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For short orbital periods and/or strong magnetic fields, angular momentum loss through magnetic 


braking can be completely suppressed. 


The reduction of magnetic braking for a given wp can be parametrized with the fraction of 
open field lines (®). In the extreme case of no wind zone, i.e. uwnp is strong enough to prevent 
any magnetic braking, ® = 0.0. On the other hand, in case the impact of the magnetic moment 
is negligible, the second dead zone does not exist and magnetic braking is not reduced, leading to 
® = 0.258. This prescription for the reduction of MB depends on the binary parameters and uwp, 
and therefore allows us to incorporate the impact of strong WD magnetic fields in binary evolution 


codes*?”°, 


The resulting angular momentum loss due to magnetic braking in the presence of a strongly 


magnetic white dwarf (JMB mag) can be written as 


, ; @ \ 
JMB,mag = JMB,non—-mag 0.258 ; (13) 


where JMB sou is the angular momentum loss due to magnetic braking in non-magnetic CVs. 


In this work we used the widely adopted RVJ prescription™, with y = 3, i.e. 


M.: RoN? (Q2\° 
—30 2 4 2 2 
JMB,non—mag = —3.8 x 10 (=) Re (=) (>) , 14) 


where Rp is the donor radius and {2 is the donor spin, which is the same as the orbital angular 
velocity. This reduced magnetic braking model has been incorporated in binary population models 
of magnetic CVs and has been shown to convincingly reproduce the orbital period distribution of 


magnetic CVs”° and we here implemented the same prescription in MESA. 
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Modelling magnetic CV evolution with MESA. The above described effects were incorpo- 
rated in the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA) (ver- 
sion r10108)458:22, The MESA equation of state is a blend of the OPAL}, SCVH*’, PTEH®!, 
HELM”, and PC® equation of states. Radiative opacities are primarily from OPAL™*:*, being the 
regime low-temperature dominated”® or high-temperature, Compton-scattering dominated’. Elec- 
tron conduction opacities are also included**. Nuclear reaction rates are a combination of rates 
from NACRE”’, JINA REACLIB®, plus additional tabulated weak reaction rates°!*?. Screening 
is included™ as well as thermal neutrino loss®°. Roche lobe radii in binary systems are computed 
using Eggleton’s fit® and mass transfer rates in Roche lobe overflowing binary systems are deter- 


mined following Ritter’s prescription®’. 


The parameters of our standard model were set as follows. During the whole CV evolution 
we assumed that all material accreted by the WD is lost through nova eruptions. In Equation 
this corresponds to setting « = 0. For the geometry of nova shells, we assumed that mass leaves 
the WD with some asymmetry (i.e. neither completely spherical nor radial), which implies values 
of 7 in Eq. [1] between 2/3 and 1. We set 7 = 2.5/3. Additional mass loss from the system®*°? 
was taken into account by using a = 0.75 in Equation [I] The donor magnetic field was assumed 
to be By = 1 kG, which is consistent with recent measurements of field strengths of low-mass 
stars”. We fixed the WD mass to 0.8 Mo, which corresponds to the average value found among 
observed CVs?*:”!, For the white dwarf magnetic field, we assume Bwp = 60 MG, which is 
slightly stronger than the average field strength of polars’, but corresponds to the average field 


strength found in pre-polars?. We assumed that the dynamo is efficient only if the core of the 
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white dwarf is crystallizing (between 10 and 80 per cent crystallized) and the spin period of the 
white dwarf is ~ 20 — 30s. Finally, for the synchronisation time-scale, we set r = 1 Myr, which 
is consistent with the expected rapid spin-down in close binaries harbouring a rapidly rotating 


magnetic white dwarf such as AR Sco?. 


Using this set of parameters our model can reproduce the observations of magnetic detached 
and semi-detached close white dwarf binary stars (Supplementary Figure 1). According to this 
model, pre-polars and AR Sco have been CVs in the past and at that time the strong magnetic 
fields of their white dwarfs were generated. This prediction is in perfect agreement with the cool 
white dwarf temperatures of these detached systems, with them being close to Roche-lobe filling, 


and with their preferred occurrence in the 3-5 hr orbital period range. 


We emphasise that no fine-tuning was required to obtain this agreement. In contrast, the 
model naturally explains the observations and variations of the key parameters do not affect the 
general predictions (Supplementary Figures 2, 3 and 4). For instance, by changing the spin-up 
parameters (i.e. a, 7) and €), the only effect is to accelerate or delay the spin-up of the white dwarf 
and slightly change the minimum spin period that can be reached”*. For all reasonable values of 
these parameters the white dwarf will manage to spin-up to sufficiently large spin rates so that the 


conditions for generating strong magnetic fields can eventually be met. 


Similarly, varying the synchronisation time-scale 7 affects the synchronising torque but does 
not change the evolutionary sequence (Supplementary Figure 2). The shorter 7, the stronger the 


spin-down torque and therefore the longer the detached phase. Only for 7 = 100 Myr, the detached 


19 


phase can be fully prevented if, in addition, the magnetic field strength of the white dwarf does not 
exceed ~ 200 MG. The latter condition arises as for very strong fields reduced magnetic braking 
alone can cause the system to detach in full analogy to the mechanism producing the classic orbital 


period gap. 


A slightly stronger constraint on the synchronisation time scale results from the fact that it 
should not exceed the duration of the generated detached phase as otherwise the observed syn- 
chronised pre-polars would remain unexplained. For white dwarfs rotating near break-up, in order 
to have a synchronised system after the detached phase, time scales cannot exceed 30 — 40 Myr. 
This maximum time scale is definitely not a problem for our model, since the synchronisation time 
scales for the range of magnetic field strengths we consider here, i.e. those of polars, pre-polars 


and AR Sco, are expected to be well below 30 — 40 Myr”. 


The overall evolution does also not sensitively depend on the assumed white dwarf and sec- 
ondary star magnetic fields. Varying these field strengths changes only slightly the condition for 
synchronisation (Equation (9). For weaker/stronger Bə, this condition becomes more/less restric- 
tive. Increasing or decreasing Bwp affects somewhat the orbital period at which synchronism 
can be obtained, which is longer for stronger fields (Supplementary Figure 3). The white dwarf 
magnetic field strengths also affects the duration of the detached phase. The stronger the white 
dwarf magnetic field, the more angular momentum loss due to magnetic braking is reduced and 


the longer it takes the system to evolve back into a semi-detached configuration. 


The largest and most critical assumption in our model certainly are the conditions we assume 
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for magnetic field generation which represent a strong but reasonable simplification. According 
to the most popular scaling law, the dynamo saturates for fast rotation (in white dwarfs for spin 
periods below 90 s?°) and the produced field strength depends then solely on the convective energy 
flux. However, as outlined above, the efficiency of the dynamo is also expected to significantly 
depend on the magnetic Prandtl number which is not considered in any currently available scaling 
law. As long as a reliable scaling law taking into account this dependency is not available, describ- 
ing the evolution of close magnetic white dwarf binaries needs to be based on simple assumptions 
such as those made in this paper. As soon as realistic scaling laws for the dynamo operating in 
crystallizing white dwarfs are available, we plan to incorporate those together with detailed white 


dwarf models into our simulations. 


The aim of this paper is to show how the evolution of close white dwarf binaries is affected if, 
and that is a reasonable assumption, crystallization and fast rotation can generate a strong magnetic 
field. We would expect the generated field strength in these models to depend on the rotation rate 
and the energy flux in the convection zone. Given the large range of orbital periods at the end of 
common envelope evolution, and the therefore large variety of evolutionary time scales towards the 
CV phase, we expect that such a more detailed model can explain the large range of field strength 
observed in magnetic CVs including those of intermediate polars. In case of the latter the field 
is simply not strong enough to synchronise the white dwarf spin with the orbital motion and the 


system does not enter a detached phase apart from the usual 2 — 3hr orbital period gap. 


In the Supplementary Information we present a detailed discussion of different evolutionary 
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channels predicted by our model and estimate the fraction of magnetic and non-magnetic white 


44,73-16 ac well as 


dwarfs in CVs taking into account important parameters for binary star evolution 
potential accretion heating and chemical impurities of the core of the white dwarf’. We also 
elaborate on the implications of our scenario for the generation of magnetic fields in hot single 
white dwarfs®*® as well as magnetic white dwarfs in wide binaries®!. We find that these magnetic 
fields are most likely generated by alternative mechanisms such as fossil fields resulting from main 


sequence star mergers®*®?, double white dwarf mergers!®, or common envelope mergers!®:**+®°, 
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Supplementary Discussion and Figures 


The fraction of magnetic white dwarfs in CVs. We have shown that a crystallization and rotation 
driven dynamo can explain the properties of all close and strongly magnetic white dwarf binaries 
and the absence of magnetic white dwarfs in young detached white dwarf binary stars. In what 
follows we discuss under which conditions our model is likely to also reproduce the observed 
fraction of strongly magnetic white dwarfs in CVs. The measured fraction of strongly magnetic 
white dwarfs in CVs is more than one third (15 out of 42 systems within 150 pc’). If the briefly 
detached pre-polars and AR Sco are included in this sample (the four pre-polars SDSS 0303+0054, 
WX LMi, SDSS J1212+0136, SDSS J1452+2045 and AR Sco are within 150 pc*!°), the fraction of 
systems with magnetic white dwarfs further increases to 20 of 47 systems corresponding to 43 + 6 


per cent. 


At first glance, our model might seem to produce too many magnetic white dwarfs in CVs as 
every white dwarf in a non-magnetic CV accretes angular momentum and is therefore expected to 
become rapidly rotating in a short period of time. In addition, if the white dwarf continuously cools 
it should start to crystallize and one might expect that our model predicts every white dwarf ina CV 
to eventually develop a strong magnetic field. However, accretion can increase the core temperature 
of white dwarfs in CVs” and prevent or delay crystallization. In addition, the generation of a strong 
magnetic field through the rotational and crystallization driven dynamo only works if fast rotation 
and a crystallizing core with an outer convection zone are present at the same time. The strength 


of the magnetic field generated by this dynamo depends on the convective energy flux” and most 
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likely also on the magnetic Prantdl number. However, in the absence of a scaling relation including 
the dependence on the latter, it remains unclear at which stages of crystallization the conditions 
are suitable for generating strong magnetic fields. In our model we assumed that in addition to fast 
rotation at least 10 per cent and not more than 80 per cent of the core need to be crystallized for 
the dynamo to be able to generate a strong magnetic field but this assumption is rather arbitrary. 
In what follows we describe the three different types of CVs that are non-magnetic according to 
our model and explain how the critical core mass ratios and other important parameters impact the 


predicted relative numbers of magnetic and non-magnetic CVs. 


First (channel 1), our model predicts the white dwarfs in CVs that did not yet accrete enough 
angular momentum remain non-magnetic. How long this phase takes, and therefore also the rel- 
ative number of CVs that are expected to be found in this stage in observed samples, depends on 
the efficiency of the accretion of angular momentum and the mass transfer rate (see equation 1 in 
the main text) which in turn depends on the angular momentum loss. For instance, CVs born with 
fully convective secondaries are expected to be driven by gravitational radiation only and typically 
have mass transfer rates 1-2 orders of magnitude lower than systems with earlier secondary stars. 
Therefore, the white dwarfs born with fully convective secondary stars will take longer to spin up 
to near-break velocities than those with earlier secondary stars that start mass transfer at longer 
orbital periods. However, as the spin up time scale is in general short compared to the secular time 
scale of CVs, we expect binary population synthesis based on our model to predict rather few CVs 


to be in this non-magnetic CV state. 
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Second (channel 2), CVs with white dwarfs that did not yet start the crystallization process in 
their cores remain non-magnetic CVs until the required condition is reached. How long this takes 
determines how many of such systems are predicted to exist in observed samples and depends on 
the fraction of crystallized core-mass that is assumed to be necessary to make the dynamo work 
(we assumed 10 per cent), on the time it takes the system to evolve from the common envelope 
to the CV phase, and whether accretion can heat up the core’”’”* to avoid or delay crystallization. 
Assuming core heating is inefficient, according to recent white dwarf models”* the time needed for 
a 0.6, 0.8 and 1.0 Mowhite dwarf to have 10 per cent of its interior crystallized is approximately 
3.2, 1.9 and 1.1 Gyr. This time span might be substantially increased for smaller white dwarf 
masses as the magnetic field needs to diffuse through a radiative zone which sits on top of the 
convective mantle. Only for massive white dwarfs the convective mantle extends close to the 
surface*’. This implies that detached binaries that emerged from common envelope evolution with 
short orbital periods, i.e. close to the start of mass transfer, can remain non-magnetic CVs for 
23 Gyr. If core heating due to accretion is efficient and prevents the white dwarf from crystallizing 
during the CV phase, these CVs will even remain non-magnetic forever. Therefore a non-negligible 
number of non-magnetic CVs descending from short period post common envelope binaries is 


predicted to exist in observed samples of CVs. 


Finally (channel 3), CVs containing white dwarfs that significantly crystallized before being 
spun up to fast rotation due to accretion, will not develop a strong magnetic field. In other words, 
if the white dwarf had enough time too cool so that a large fraction of its core is crystallized before 


the secondary fills its Roche-lobe, the convective energy flux in the outer core will not be sufficient 
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to generate a strong magnetic field. These systems will remain non-magnetic forever under the 
assumption that core heating due to accretion is inefficient. Even if accretion heating can melt a 
significant fraction of the core, this process is supposed to take several Gyr”, i.e. the system would 
remain non-magnetic for a large fraction of its CV lifetime. Thus, whether a given system becomes 
a non-magnetic CV formed through this channel depends mostly on the time it spends between the 
white dwarf formation, i.e. the common envelope phase, and the onset of accretion. How many 
systems of a predicted CV population spend enough time in this detached white dwarf binary 
phase depends sensitively on the crystallized mass fraction (we assumed 80 per cent) at which 
the dynamo becomes unable to generate a strong magnetic field, on the inter-correlated initial 
orbital period, mass ratio, and eccentricity distributions of the progenitor systems, on the common 


envelope efficiency, and on the strength of angular momentum loss due to magnetic braking. 


While the relative number of non-magnetic CVs predicted by our model through channel 1 
is most likely very small, we expect that the other two groups of systems can easily explain the 
~ 50 — 70 per cent of non-magnetic systems in the observed CV sample. Which of the two 
scenarios contributes more to the fraction of non-magnetic CVs depends on several assumptions. 
For example, for a large common envelope efficiency, relatively inefficient magnetic braking, and 
a critical crystallized mass fraction below 80 per cent, the majority of systems will evolve through 
channel 3 and will remain non-magnetic for at least several Gyr. In contrast, for a small common 
envelope efficiency, efficient magnetic braking, a large value for the critical fraction of crystallized 
core mass, and inefficient core heating, most CVs will eventually become magnetic but many of 


them may spend several Gyr as non-magnetic systems before the dynamo starts operating (channel 
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2). 


The parameters (strengths of magnetic braking, critical crystallized mass fractions, common 
envelope efficiency, efficiency of core heating) are essential for predicting the evolution of close 
white dwarf binaries and for determining the fraction of strongly magnetic CVs predicted by our 
model. Unfortunately, several of these parameters are currently neither theoretically nor obser- 
vationally well constrained. The common envelope efficiency is most likely relatively small as 
indicated by observational surveys of detached white dwarf binaries’? and heating of the core is 


77.78 In contrast, the strength of magnetic braking is unknown® and 


likely to play an important role 
the same is true for the critical mass fractions (as detailed simulations for the dynamo we propose 
here do not yet exist). In addition, the initial binary parameter distributions usually assumed in CV 
population models are uncorrelated which is inconsistent with observational surveys performed 
in the past decades™®. These surveys show that the orbital period, primary mass, mass ratio and 
eccentricity distributions are inter-correlated. This is important because the combination of these 
parameters determines whether a given binary system will survive or not the common-envelope 
phase and which of those that survive and become CVs within the Hubble time, have relatively 


long or short detached lifetimes (possibly becoming slightly or highly crystallized) before the on- 


set of mass transfer. 


Based on binary population models we performed in the past’*’’°, using standard parameters 


for CV evolution (common envelope efficiency of 0.25, magnetic braking according to Rappaport 


1.“4 


et al.“*, assuming the critical crystallized mass fractions used in this paper (i.e. 10 and 80 per cent), 


28 


and ignoring core heating of the white dwarf through accretion, we estimate the predicted fraction 
of magnetic systems to be roughly 50 per cent. While this value slightly exceeds the currently ob- 
served number of magnetic CVs, relatively small changes in the prescription for magnetic braking 
and/or the critical mass fractions of crystallized matter in the core of the white dwarf can easily 
bring into agreement the observed fractions and the predictions of our model. For the same as- 
sumptions but including efficient core heating, the predicted fraction of magnetic CVs decreases to 
~ 20 per cent. This estimate is below the measured fraction but our conclusion remains unchanged 
as relatively small changes in uncertain parameters such as critical core mass and magnetic braking 


could increase the number of magnetic systems to match the observed fraction. 


For the sake of completeness, as an additional unknown parameter, we emphasize that also 
the metallicity of the white dwarf progenitor star could play a role since the settling of neutronized 
impurities like 22Ne and 56Fe caused by solidification can release more energy than the separation 
of the C/O mixture which significantly affects the cooling of the white dwarf and the crystallization 


process”. 


As a final remark, we state that the small number of observationally discovered pre-polars 
and the even smaller number of AR Sco-like systems are in agreement with the relatively short 


time scale that our model predicts for these systems. 


Implications for single white dwarfs and white dwarfs in wide binaries. The model presented 
in this paper describes a new evolutionary sequence for close white dwarf binaries which explains 


the origin and the relative numbers of strongly magnetic white dwarfs in these systems. While 
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solving the puzzle for the occurrence of strongly magnetic white dwarfs in close binaries represents 
a key result on its own, we here briefly mention the implications for single white dwarfs and white 


dwarfs in wide binary stars. 


It is clear that the dynamo proposed here to operate in CVs cannot explain all strongly mag- 
netic single white dwarfs because it requires fast rotation and crystallization. A non-negligible 
number of strongly magnetic single white dwarfs” are hot and therefore most likely not in the pro- 
cess of crystallization and, even worse, several of the highest field single white dwarfs are known 
to have rather long spin periods’, 80. While field strengths of the order of 0.1 — 1 MG might be 
generated in crystallizing single white dwarfs with rotation periods of the order of hours or days, 
strong magnetic fields in single white dwarfs must be generated at least partly by a different mech- 
anism. The same conclusion holds for the recently discovered strongly magnetic white dwarfs in 
wide binaries®'. In most cases the separation between the two stars in these binaries exceeds sev- 
eral hundred astronomical units and therefore a previous phase of mass transfer and spin-up of the 


white dwarf can be excluded. 


Similarly, the main mechanism for generating these strong magnetic fields in young single 
white dwarfs and wide binary stars must not be operating in close binaries as we would other- 
wise observe strongly magnetic white dwarfs in young detached close white dwarf binaries. The 
candidate scenarios that may fulfill this conditions are white dwarf mergers!’, fossil fields!’, or 
the common envelope merger scenario'®*+°°, Given that the most likely progenitors of magnetic 


white dwarfs according to the fossil field scenario are main sequence star mergers®”"®?, all three 


30 


scenarios require single magnetic white dwarfs to form from binary stars. This implies that all 
scenarios require initial triple star configurations to produce strongly magnetic white dwarfs in 
binaries. As triples are intrinsically less common than binary stars, all three models may produce 
more strongly magnetic single white dwarfs than strongly magnetic white dwarfs in close and wide 
binaries. Which of the mechanisms is most suitable to explain magnetic single white dwarfs and 
magnetic white dwarfs in wide binaries depends on the relatively poorly constraint distribution of 


initial triple star configurations. 


We conclude that the crystallization and rotation driven dynamo proposed here explains the 
large number of magnetic white dwarfs in CVs, as well as the existence and characteristics of 
AR Sco and the pre-polars, and might be complemented by e.g. the double white dwarf merger 
and/or the common envelope merger scenario which are good candidates for explaining strongly 


magnetic single white dwarfs. 
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Supplementary Figure 1. Comparison of different evolutionary tracks with observations. The de- 
tached phases (dashed lines) explain the existence of only cool magnetic white dwarfs in detached 
binaries that are close to Roche-lobe filling because the magnetic field was generated during the 
CV phase when the white dwarf was old and crystallizing. The model also offers an explanation 
for AR Sco which in fact is a progenitor of the pre-polars. In AR Sco, nicely reproduced by the 
blue track, the secondary detached from its Roche-lobe but spin and orbital motion are not yet 
synchronised. The subsequent detached synchronised phase explains the orbital periods, Roche- 
lobe filling factors, and white dwarf temperatures of pre-polars. To generate the tracks we only 
changed the initial orbital period (P,) and initial secondary mass (Mə), as well as the white dwarf 


spin period (Pwp) when the magnetic field is generated. The parameters we used were the fol- 
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lowing P, = 0.14days, Mz = 0.17Mo, Pwo = 30s (purple), P, = 0.68days, My = 0.55Mo, 
Pwo = 15s (blue), P, = 0.925 days, Mz = 0.75Mo, Pwo = 20s (green), P, = 1.544 days, 
M> = 1.2Mo5, Pwo = 30s (orange). In all tracks the white dwarf mass was assumed to be 
0.8M.”° and the generated magnetic field strengths was assumed to be 60 MG. The errorbars on 


the observations correspond to lo. 
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Supplementary Figure 2. Evolutionary sequences highlighting the impact of different choices of 
the synchronisation time scale. The shorter the time scale, the more pronounced and the longer 


the detached phase. For the longest time scales (blue, purple), accretion starts before the system 


synchronised. Time scales exceeding 1 Myr are most likely not realistic. 
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Supplementary Figure 3. Evolutionary sequences highlighting the impact of the white dwarf 
magnetic field strength. The stronger the generated field, the longer the detached phase as angular 
momentum loss due to magnetic braking is more reduced. For the weakest field strengths (brown, 


red), a less pronounced standard gap still occurs when the secondary becomes fully convective. 
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Supplementary Figure 4. Evolutionary sequences highlighting the impact of different choices of 
the white dwarf spin period (Pwp) at onset of magnetic field generation. The hatched region indi- 
cates the locus of the standard orbital period gap. In each track, solid and dashed lines correspond 
to semi-detached and detached phases respectively. With the exception of the case Pwp = 14s 
(brown curve), a detached phase is generated above the standard orbital period gap, at an orbital 
period ~ 3.75hr. For Pwp at least twice the break-up period, magnetic field generation starts 
at the same period and the slower the spin, the shorter the detached phase. To spin-up the white 
dwarf to velocities close to the break-up period (brown track), the spin-up phase is much longer 


and therefore the pronounced detached phase occurs at a significantly shorter orbital period. 
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